Priedel oscillations and the Kondo screening cloud 



00 
O 

o 

Oh: 

< 

o 



CO 

c3 



i 

c 

o 
o 



> 
o 

oo 

(N 

o 

(N 
O 
OO 
O 



Ian Affleck 1 , Laszlo Borda 2 and Hubert Saleur 3,4 
1 Department of Physics and Astronomy, University of British Columbia, Vancouver, B.C., Canada, V6T 1Z1 
2 Physikalisches Institut, Universitdt Bonn, Nussallee 12, D-53115 Bonn, Germany 
3 Service de Physique Theorique, CEA Saclay, Gif Sur Yvette, 91191, France 
4 Department of Physics and Astronomy, University of Southern California, Los Angeles, CA 90089, USA 

(Dated: April 30, 2008) 

We show that the long distance charge density oscillations in a metal induced by a weakly coupled 
spin-1/2 magnetic impurity exhibiting the Kondo effect are given, at zero temperature, by a universal 
function F(r/£ic) where r is the distance from the impurity and fx, the Kondo screening cloud size 
= Kvf /(fcsTif), where vf is the Fermi velocity and Tk is the Kondo temperature. F is given by a 
Fourier-like transform of the T-matrix. Analytic expressions for F(r /£k) are derived in both limits 
r <C fx and r ^> £k and F is calculated for all r/£zc using numerical methods. 



The interaction of a single magnetic impurity with the 
conduction electrons in a metal is often described by the 
Kondo model: 



H = J2 e(k)4j, a + J £ JlJfifof, • S. (1) 
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In general we also include a potential scattering term 
in the Hamiltonian: H^H + V^Tzz, W- ipz, ■ This 
model exhibits a remarkable crossover from weak to 
strong coupling behavior as the energy scale is lowered 
through the Kondo temperature, ksTx ~ 2?exp(— 1/Ao) 
where V is an ultra-violet cut-off scale (such as a band 
width) and Ao is the dimensionless bare coupling constant 
(= Jv where v is the density of states, per spin). For a re- 
view, see, for example, chapter 4 of |l| . The renormalized 
Kondo coupling, A(E), becomes of 0(1) at E ~ UbTk- 
While physics at energy scales E 3> ^bTk is given by 
weak coupling perturbation theory, at E <C fesTV, the 
physics is governed by the strong coupling fixed point 
corresponding to a screened impurity and a tt/2 phase 
shift for the low energy quasi-particles. 

The length dependence of Kondo physics is much less 
well understood. It is generally expected that physical 
quantities exhibit a crossover at a length scale = 
hvp I '{ksTx) (where vp is the Fermi velocity), which is 
typically in the range of .1 to 1 micron. (We henceforth 
set h and ks to 1.) For a review see section 9.6 of [l|. 
See @, i H for original work on the subject. However, 
such a crossover at this long length scale has never been 
observed experimentally and has sometimes been ques- 
tioned theoretically [5| . One way of observing this length 
scale is through the density oscillations around a mag- 
netic impurity 0, 0, @, 0] ■ It was pointed out in that 
these should only approach the standard Friedel form at 
distances r » (k, with a form at shorter distances con- 
trolled by the T-matrix. However, experimental data so 
far doesn't seem to support this expectation 0], yield- 
ing much shorter characteristic lengths. (See also @, @].) 
One purpose here is to present a more complete theoreti- 
cal treatment of these density oscillations, since scanning 
tunnelling microscopy of magnetic ions on metallic sur- 



faces provides a new experimental technique by which 
they might now be measured. Alternative approaches 
to observing this fundamental length scale involve ex- 
periments on mesoscopic structures with dimensions of 

0((k) |. 

We focus on the case of an S = 1/2 impurity, and 
a spherically symmetric dispersion relation (normally 
e(k) = k 2 /2m— ef)- We consider this model in dimension 
D = 1, 2 or 3. 

There are two reasons why one might be skeptical that 
the length scale £k would show up in the charge density. 
One is the idea of "spin-charge" separation in D = 1 . 
The Hamiltonian of Eq. |lj in any dimension, can be 
mapped into a ID model by expanding in spherical har- 
monics and using the fact that only the s-wave harmonic 
interacts with the impurity in the case of a (5-function 
interaction. The low energy degrees of freedom of non- 
interacting ID electrons can be separated into decoupled 
spin and charge excitations, using bosonization. It is 
possible to write the Kondo interaction in terms of the 
spin degrees of freedom only and hence, one might ex- 
pect the charge density to be unaffected by the Kondo 
interaction. The fallacy in this argument is that the 
charge density at location r in the ID model contains 
a term ipL a ( r )' t pRa( r ) exp(— 2ikpr) + h.c. where R and L 
label right and left movers. Standard bosonization meth- 
ods imply that this term involves both spin and charge 
bosons: sin( \f2~n(j) c + 2kF r ) cos(y/2n(j) s (r)). This is unlike 
the term V4, a (O^Atc* ( r ) which only involves the charge 
boson. 

Another reason why one might expect no interest- 
ing Friedel oscillations follows from consideration of the 
particle-hole (p-h) symmetric case. This symmetry is 
exact, for example, in a nearest neighbor tight-binding 
model at 1/2-filling with the Kondo coupling occuring 
at the origin only. Then it can easily be proven that 
< ipj a i/jja >= 1 f° r all s it cs 3- However, a realistic 
model always breaks particle-hole symmetry. This can 
be achieved by taking a non p-h symmetric dispersion 
relation - for instance moving the density away from 1/2- 
filling in the tight-binding model. Alternatively, poten- 
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tial scattering can be included in the model. Then, p-h 
symmetry is broken even if the dispersion relation doesn't 
break it. 

We find for the density oscillations at zero temperature 
and r ^> l/k F : 



P(r)-f>a -> ^§[cos(2k F r-nD/2 + 26 P )F(r/^ K ) 



cos(2k F r - irD/2)]. 



(2) 



Here F(t/£k) is a universal scaling function which is the 
same for all D, Sp is the phase shift at the Fermi sur- 
face produced by the potential scattering, G3 = l/(47r 2 ), 
C2 = l/(27r 2 ) and C\ = 1/(2%). In general, there are 
non-zero oscillations but they vanish exactly in the p-h 
symmetry case for D=l where dp = 0, k F = it/2 and r 
is restricted to integer values, corresponding to a tight- 
binding model at 1/2-filling. In the limit of zero Kondo 
coupling, F = 1 and we recover the standard formula for 
Friedel oscillations produced by a potential scatterer (in 
the s-wave channel only). For a small bare Kondo cou- 
pling, Ao <C 1, F(t/^k) is close to 1 at r <C £,k so that the 
oscillations are just determined by the potential scatter- 
ing, oc cos(2fci?r — nD/2 + 28p) — cos(2kpr — irD/2), van- 
ishing if Sp is also zero. However, at r 3> £,k, we find that 
F(r/^K) — > —1 which is equivalent to Sp — > Sp + n/2. 
We again recover the potential scattering result but now 
the phase shift picks up an additional contribution of 7r/2 
from the Kondo scattering. 

To derive these results, following it is convenient 
to relate the scaling function, F(r/£x) to the T-matrix, 
T(w) which has already been well-studied by a number 
of methods and is a universal scaling function of lo/Tk- 
This can be done using the standard formula for the (re- 
tarded) electron Green's function: 

G(?J,w) = G {r-r<,u) + G {r,u)T(u>)G {-r",u), 

(3) 

where Go is the Green's function for the non-interacting 
case (with J = V = 0). This result is a direct conse- 
quence of the assumed S- function form of the Kondo (and 
potential scattering) interaction. The density is obtained 
from the retarded Green's function by: 



2 f° 

p(r) = / dui Im G(r,r,uj). 

t J -00 



(4) 



(The factor of 2 arises from summing over spin.) The 
exact non-interacting Green's function is: 



Go — 



—ikp 
vpk 



-I {D-l)/2 

-ik 



2nr 



exp[ifcr], (D = 1, 3) 



= -[k F /{nv F )]K -ikr , (D = 2) 



(5) 



where k = \Jkp + 2k F Lu/v F and K (z) is the modified 
Bessel function. This gives the asymptotic behavior at 



r > l/k F , lu-^V: 
Gl(r,u) -> — l - 



-iki 



n D-l 



2irr 



exp[2ik F r+2iLur/v F ]. (6) 



(This asymptotic behavior holds for general dispersion 
relations.) The T matrix in D-dimensions can be writ- 
ten at ui <C T> : T(u) = t(oj /Tk) / (2ttvd) where t is a 
universal dimensionless function of uj/Tk and vpn the 
density of states per spin at the Fermi energy, has the 
value vp) = kp" 1 /(cpiVp), with C3 = 2tt 2 , C2 — 2ir and 
C\ = 7r. Note that GqFGo is proportional to the dif- 
ference between the s-wave Green's function with and 
without the Kondo and potential scattering interactions, 
since the other spherical harmonics are unaffected by the 
interactions and cancel in G — Go- The effect of the 
s-wave potential scattering at long distances is just to 
multiply the s-wave Green's function by the phase e 2lSp , 
thus giving: 



t(u/T K ) = e 2t5p [t K {uo/T K ) 



(7) 



where tKioj /Tk) is the part of the t-matrix coming from 
the Kondo scattering. Combining Eqs. |3][7|) gives: 



p{r) 



Po 



{c D /[TT 2 v F (2Trr) D - 1 ]}Im 




D-l e 2ik F r 



x / dwe 2wr l VF [{t K {io/T K ) + i)i 



2iS r 



.(8) 



Essentially this formula (for D = 3 only) was derived 
in Q except that our treatment of p-h symmetry break- 
ing is quite different. Furthermore we apply much more 
complete knowledge of the T-matrix. We expect this for- 
mula to be valid whenever r 3> l/k F , regardless of 
the ratio r/^K- The function tuiy /Tk) is determined 
from the p-h symmetric Kondo interaction and so obeys: 
t* K (iO /Tk) = — tic(— w/Tk). Furthermore tiui/Tx) is an- 
alytic in the upper half complex uj plane since it is ob- 
tained from the retarded Green's function. It then fol- 
lows that J_ duo eKp(2iujr/v F )tK(ijj/Tji) is purely real. 
A rescaling of the integration variable implies that we 
may write: 

due 2iur ^t K (oj/T K ) = [v F /(2r)][F(rT K /v F ) - 1], 

(9) 

where the universal scaling function F is purely real. 
Thus: 

p (r)-p ^{c D /[2n 2 (2Tr) D - 1 r D ]} 

xlm {(-i) D -^ e 2ikFT [F(r/U)e 2lSp - l] } , (10) 

giving the result announced in Eq. (|2|) . While this deriva- 
tion assumed that the Kondo interaction is a spatial 5- 
function leading to the simple result, Eq. ([3]), we expect 
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our asymptotic formula for p(r) to be much more gener- 
ally true, at length scales large compared to the range of 
the Kondo interaction. 

A perturbative calculation of the T-matrix [lCj gives: 



t K (uj) = -(3i7r 2 /8)[A 2 + Aq ln(2?/w) 2 + 



fill 



where T> is of order the ultraviolet cut-off. The quan- 
tity in brackets can be recognized as the first 2 terms 
in the expansion of the square of the running coupling 
A 2 (w). For uj > T K , A(V) -> 1/ h\(\u\/T K ), so one 
expects tx — > — 37r 2 i/[8 ln 2 (|u;|/Tft-)]. Substituting the 
perturbative expansion into Eq. (J9]), gives: F{t/^k) = 
1 - (3tt 2 /8)[Ao + 2A^ln(r/a) + ...], where a is a short 
distance cut-off of order vpjT>. Again, we recognize the 
first terms in the expansion of A 2 (r), implying the short 
distance behavior: 

W60->l-3^ 2 /[81n 2 (Wr)], (r « 60- (12) 

It is an interesting fact that F{v/^k) is apparently 
given by renormalization group improved perturbation 
theory for r <C £k ■ This is quite unlike the situation for a 
related quantity, the Knight shift Q • This is again given 
by a scaling function, x( r ) = (l/?A')/( r /C-R") a t zer0 tem- 
perature. However, in this case the term of O(Ag) has a 
coefficient which diverges as the temperature T — ► 0, 
even at a fixed small r. This means that the Knight shift 
at short distances (r -C £_Ff) is not given by renormal- 
ization group improved perturbation theory, unlike the 
Friedel oscillations. Instead, the Knight shift exhibits a 
non-perturbative behavior, even at short distances. A 
conjecture was made for this non-trivial short distance 
behavior in 0]. The fact that F(r/^K) is perturbative at 
small r/^K seems to follow from the fact that T(lo/Tk) is 
perturbative at large lo/Tk together with Eq. (J9j> , which 
presumably implies that the short distance behavior of F 
is given by the high-frequency behavior of T(uj/Tk)- The 
general question of which quantities are perturbative or 
non-perturbative at short distances in the Kondo model 
(and other quantum impurity models) remains open. 

Perturbation theory for the Friedel oscillations breaks 
down at r of 0(£k) but at r 3> £,k we may use Nozieres 
local Fermi liquid theory. This gives the T-matrix: tx — > 
-i[2 + iuj/T B - 3ui 2 /4T B + ...]. Here T B corresponds to 
a particular definition of the Kondo temperature. (See, 
for example, chapter 4 of |l].) It is related to the Wilson 
definition, called simply Tk in l] by Tp — 2Tk/(ttw) 
with the Wilson number w s» .4128. Substituting in Eq. 
© gives: 



(r » 6f ), 



3(™m/(32r 2 



(13) 



where we have defined £k precisely in terms of the Wilson 
definition of Tk- £,k = ^f/Tk- Nozieres' perturbation 
theory can be turned into a full perturbation theory 11 1 
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FIG. 1: NRG results on charge oscillations around a Kondo 
impurity coupled to ID conduction electrons with particle- 
hole symmetry. Note that the oscillations vanish at kpr/n g 
N. As shown in the inset, the properly rescaled envelope 
function of the oscillations (extracted as p — pa at the local 
maxima) for different Kondo couplings nicely collapse into 



one universal curve except for the points where r ~ k. 



In 



the inset we show the analytical results for the asymptotics as 
well: Note the good agreement between the analytical results 
and the numerics. 



by taking into account more irrelevant operators in the 
vicinity of the low energy fixed point, which give higher 
order terms in Eq. (|13p . 

In order to strengthen our analytical results, we have 
performed extensive numerical renormalization group 



(NRG) calculations [12j,ll3|. In Wilson's NRG technique 
-after a logarithmic discretization of the conduction elec- 
tron band- one maps the original Kondo Hamiltonian to 
a semi-infinite chain with the impurity at the end. As a 
direct consequence of the logarithmic discretization the 
hopping amplitude along the chain falls off exponentially. 
This separation of energy scales allows us to diagonalize 
the chain Hamiltonian iteratively in order to approxi- 
mate the ground state and the excitation spectrum of 
the full chain. If one is interested in spatial correlations, 
however, some care is needed. The cornerstone of the 
model, the logarithmic discretization causes not only the 
exponential fall-off of the hopping amplitude, but also a 
very poor spatial resolution away from the impurity. To 
tackle that problem, we introduce Wannier states cen- 
tered both around the impurity and the point of interest 
r thus reducing the problem to a two impurity type calcu- 
lation. Such an approach has been demonstrated to work 
recently by evaluating the spin-spin correlation function 
around a Kondo impurity -see 14| . To get the amplitude 
of the charge oscillations one needs the explicit value of 
Uf which we obtained by calibrating the NRG code with 
a pure potential scattering model. 

We show results for different Kondo couplings in Fig. 
flT]). p(r) — p Q ~ sm(2kpr) in agreement with Eq. @ 
for Sp = 0, the expected p-h symmetric result since 
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we use a flat symmetric band with no potential scat- 
tering. In the inset of FigQ] we show NRG results for 
F{t/^k) showing good agreement with the asymptotic 
predictions of Eqs. (|12[) and (fT5|) and fair agreement with 
the prediction of the "one spinon approximation" 0, UH j 
t K = -2i/[l - iuj/T B ], F(u) = 1 + Auae 2ua Ei (-ua), 
a = T B /T K = 2/(vrw) sa 1.542. (Ei is the exponential- 
integral function.) This is a challenging NRG calcula- 
tion since universal behavior is only expected to occur 
for £k , r 3> kp 1 (i.e. at distances beyond several peri- 
ods of the density oscillation and at weak coupling). On 
the other hand, the numerical error increases at large r. 
The non-universal, coupling dependent part of the charge 
density oscillations is much more extended in space than 
that of the spin-spin correlator computed in [l4|. That 
is the main source of the scattering of data points in 
the inset of Fig. (fT]). It is interesting to note, from the 
figure, that F rs 0, corresponding to the midpoint of 
the crossover from weak to strong coupling, occurs at 
r pa (0.12 ± 0.02)£/f . Thus an experimental detection 
of the Kondo screening cloud via the density oscillations 
would "only" need to measure out to distances of order 

/TO to see at least half of the crossover. In STM exper- 
iments the most readily accessible measure of the Kondo 
temperature is the half- width of Im T(cu), T\i<i sa 2Tk 
[HI]. Once this number is determined experimentally, 
then the midpoint of the cross-over of the Friedel os- 
cillations is predicted to occur at r sa vf/[5Ti/ 2 ]- At fi- 
nite temperature, Friedel oscillations decay exponentially 
with a thermal correlation length £t = 2itvf/T so it is 
necessary to be at sufficiently low T that > £j° to mea- 
sure the Kondo screening cloud. Direct electron-electron 
interactions, ignored in the Kondo model, can also lead 
to decay of the Friedel oscillations with a decay length 
related to the inelastic scattering length. However, Fermi 
liquid theory (typically believed to be valid in D=2 or 3) 
implies that this length also diverges as T — > 0. 

The Kondo screening cloud does not show up in the 
energy resolved density of states, — (2/w) Im G(r,r,u>), 
measured in STM and given by Eq. ([3]). This has a 
trivial r-dependence l/r- 0-1 at r > 1/kp. At fixed r 
the Kondo scale only enters through the w-dependence. 
Only after doing the w-integral, to get the total electron 
density does the Kondo scale appear in the r-dependence. 

Previous attempts Q to fit experimental data on den- 
sity oscillations around Cu and Mn impurities in Al to 
formulas like Eq. ([2|) have yielded characteristic lengths 
which are much smaller than £r- as determined from the 
experimentally measured Kondo temperature. We think 
these issues deserve revisiting, using STM. Im T, mea- 
sured from the energy-resolved density of states (at a 
fixed location near the impurity) , has a peak with a width 
identified as Tk- This identification is not completely 
obvious since it is typically not feasible to raise the tem- 
perature past Tk (due to diffusion of the impurity) nor 
to apply magnetic fields corresponding to Zeeman ener- 



gies of 0{Tk)- It follows from Eq. ((4]) that there should 
be a change in the envelope of the density oscillations at 
the corresponding length scale vf/Tk- An accurate mea- 
surement of p(r), if it agrees with our results, would both 
resolve an open fundamental question in Kondo physics 
and firmly establish that these systems really do exhibit 
the Kondo effect. We emphasize that the large size of 
the Kondo cloud makes it very hard to observe. At such 
large distances that F(r/£x) has changed significantly 
from its short distance asymptote of one, the l/r D fac- 
tor in Eq. @ makes the oscillations very small. Clearly 
the situation is improved in 2 dimensional systems. 

In conclusion, we have shown that the Friedel oscilla- 
tions around a Kondo impurity exhibit a universal be- 
havior characterized by the length scale We have 
determined the corresponding universal scaling function 
analytically in both limits r <C S,k and r ^> £k and 
numerically at intermediate r/^K- It exhibits renormal- 
ization group improved weak coupling behavior at short 
distances, quite unlike the Knight shift, raising intriguing 
general questions about which quantities are perturbative 
and which are not in this limit, for this and other models. 
The envelope of the oscillations, given in Eq. ([2|), exhibits 
a crossover from short to long distances corresponding to 
an increase of the s-wave phase shift by tt/2. However, 
at intermediate distances, the result does not correspond 
to simple potential scattering for any value of the phase 
shift. We have determined precisely the distance at which 
the crossover occurs in terms of the measure of the Kondo 
temperature accessible to STM experiments. 
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